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In  this  report  the  development  of  a  new  model  for  simulating  flood  inundation  is 
briefly  outlined  (see  table  1).  The  model  is  designed  to  operate  with  high  resolution 
raster  Digital  Elevation  Models  which  are  becoming  increasingly  available  for  many 
lowland  floodplain  rivers  and  is  based  on  what  we  hypothesise  to  be  the  simplest 
possible  process  representation  capable  of  simulating  dynamic  flood  inundation.  This 
consists  of  a  one-dimensional  kinematic  wave  approximation  for  channel  flow  solved 
using  an  explicit  finite  difference  scheme  and  a  two-dimensional  diffusion  wave 
representation  of  floodplain  flow.  The  model  is  applied  to  a  35  km  reach  River 
Meuse  in  the  Netherlands  using  only  published  data  sources  (see  Table  2)  and  used  to  . 
simulate  a  large  flood  event  which  occurred  in  January  1995.  This  event  was  chosen 
as  air  photo  and  Synthetic  Aperture  Radar  (SAR)  data  for  flood  inundation  extent  are 
available  to  enable  rigorous  validation  of  the  developed  model  (see  Figure  1).  100,  50 
and  25  m  resolution  models  were  constructed  and  compared  to  two  other  inundation 
prediction  techniques:  a  planar  approximation  to  the  free  surface  and  a  relatively 
coarse  resolution  two-dimensional  finite  element  scheme.  The  model  developed  in 
this  paper  out-performs  both  the  simpler  and  more  complex  process  representations, 
with  the  best  fit  simulation  correctly  predicting  81.9  %  of  inundated  and  non- 
inundated  areas.  This  compares  to  69.5  %  for  the  best  fit  planar  surface  and  63.8  % 
for  the  best  fit  finite  element  code.  However,  when  applied  solely  to  the  7  km  of  river 
below  the  upstream  gauging  station  at  Borgharen  the  planar  model  performs  almost  as 
well  (83.7  %  correct)  as  the  raster  model  (85.5  %  correct).  This  is  due  to  the 
proximity  of  the  gauge  which  acts  as  a  control  point  for  construction  of  the  planar 
surface  and  the  fact  that  here  low  lying  areas  of  the  floodplain  are  hydraulically 
connected  to  the  channel.  Importantly  though  it  is  impossible  to  generalise  such 
application  rules  and  thus  we  cannot  specify  a  priori  where  the  planar  approximation 
will  work.  Simulations  also  indicate  that,  for  this  event  at  least,  dynamic  effects  are 
relatively  unimportant  for  prediction  of  peak  inundation.  Lastly,  consideration  of 
errors  in  typically  available  gauging  station  and  inundation  extent  data  shows  the 
raster  based  model  to  be  close  to  the  current  prediction  limit  for  this  class  of  problem. 


Figure  1.  Time  series  of  inundation  extent  predicted  by  a  dynamic  simulation  of  the 
raster  model  using  a  25  m  resolution  DEM  for  the  7  km  reach  downstream  of  the 
Borgharen  gauging  station.  This  is  compared  to  the  air  photo  derived  shoreline 
sampled  at  approximately  160  hours  into  the  simulation. 
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Data  requirement 

Source 

Comments 

Raster  Digital  Elevation  Model. 

Typically  derived  from  air 

photogrammetry  or  airborne  laser 
altimetry  (LiDAR). 

Grid  resolutions  of  approximately  25- 
100m  would  seem  appropriate  for  most 
floodplain  applications,  although  smaller 
resolution  are  obviously  preferable. 
Vertical  accuracy  of  the  DEM  should 
generally  be  less  than  ±0.25  m. 

Inflow  discharge  hydrograph. 

Gauging  station  records.  Flow  enters 
the  model  through  the  upstream  channel 
cell  forming  the  first  location  on  the 
local  drainage  direction  map. 

Model  can  be  used  in  either  steady  state 
or  dynamic  modes,  but  flows  should  be 
accurate  to  ±10  %.  For  dynamic 
simulations,  temporal  resolution 

depends  on  the  speed  of  the  hydrograph 
rise  but  typically  at  least  hourly  data  are 
required. 

Channel  slope. 

Taken  from  the  DEM  or  surveyed  cross 
sections. 

Can  be  set  individually  for  each  grid  cell 
if  necessary. 

Channel  width. 

Taken  from  the  DEM  or  surveyed  cross 
sections. 

Can  be  set  individually  for  each  grid  cell 
if  necessary.  Need  not  be  the  same  as 
the  model  grid  resolution. 

Bankfiill  depth. 

Taken  from  the  DEM  or  surveyed  cross 
sections. 

Can  be  set  individually  for  each  grid  cell 

if  necessary. 

Initial  estimate  of  channel  flow  depth 

Reasonable  value  based  on  experience 
and  exmination  of  surveyed  cross 
sections  and  rating  curves. 

Model  is  run  with  constant  in-bank 
discharge  for  a  start  up  period  to  allow 
realistic  channel  water  depths  and  flow 
velocities  to  develop.  Start  up  period 
should  be  based  on  the  time  taken  for 
flood  waves  to  cross  the  domain. 

Channel  and  floodplain  friction. 

User  defined  parameters  typically 
chosen  with  reference  to  published 
tables  such  as  those  given  by  Chow 
(1959)  or  Acrement  and  Schneider 
(1984). 

Nc  typically  between  0.01  and  0.04 

Nfp  typically  between  0.03  and  0.15 

Can  be  set  individually  for  each  grid  cell 
if  necessary. 

Model  time  step 

User  defined  based  on  Courant  number 
stability  constraints.  An  explicit 
numerical  scheme  is  used  so  the  stability 
is  a  function  of  the  cell  dimensions  and 
the  flow  rate.  As  water  enters  the  model 
via  a  single  inflow  cell  at  the  head  of  the 
reach,  flow  rates  in  this  cell  are  usually 
the  limiting  factor. 

Varies  between  applications  but  typical 
values  are  in  the  range  2-20  s. 

Table  2:  Summary  of  model  data  and  parameter  requirements. 
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